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Chapter 1 

Pairing fluctuations and gauge symmetry restoration in 
rotating superfluid nuclei 

Yoshifumi R. Shimizu 

Department of Physics, Faculty of Sciences, Kyushu University, 
Fukuoka 812-8581, Japan 

Rapidly rotating nuclei provide us good testing grounds to study the 
pairing correlations; in fact, the transition from the superfluid to the 
normal phase is realized at high-spin states. The role played by the 
pairing correlations is quite different in these two phases: The static 
(BCS like mean-field) contribution is dominant in the superfluid phase, 
while the dynamic fluctuations beyond the mean-field approximation are 
important in the normal phase. The influence of the pairing fluctuations 
on the high-spin rotational spectra and moments of inertia is discussed. 

1. Rotation and pairing correlations — Introduction 

In all the different contributions in this Volume, various aspects of the 
pairing correlations, which play important roles in nuclear physics, are dis- 
cussed. In this contribution I would like to concentrate on the effect of the 
pairing fluctuations in rapidly rotating nuclei,^ which is generic and yet far 
from trivial. Here the pairing fluctuations mean the dynamic motions of the 
pairing gap, i.e., so-called the pairing vibrational (see also Ref. |3]), whose 
effects appear beyond the static (BCS) mean-field approximation and are 
characteristic in the finite system like atomic nucleus. 

It is well known that most of non-closed shell nuclei, which have 
quadrupole deformed shape, exhibit collective nuclear rotations.^ In the 
80's, the combined developments of the heavy-ion accelerators and the 
high-resolution 7-ray detectors made it possible to explore the properties of 
rapidly rotating nuclei, i.e., the high-spin states up to spin values / « 60ft, 
of medium and heavy nuclei. Many interesting phenomena and issues have 
been revealed; see e.g. Refs. [5HH and Refs. 18191 for more recent progress. 
The pairing correlations, either static or dynamic, play a crucial role in 
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most of these phenomena over the wide spin range. As for a well-known 
example, the superfluidity is responsible for the reduction of the moment of 
inertia for the collective nuclear rotation near the ground statepS it takes 
only about (or even less than) half of the rigid-body value, which is ex- 
pected for the independent nucleonic motions in the deformed mean-field. 
In this way, the ground states of deformed nuclei can be well described by 
the BCS theory with finite pairing gaps, A ss 1 MeV, and the BCS quasi- 
particles appear as a basic excitation mode. In fact, the "backbending" 
phenomenon,^ which is systematically observed at spin / « 10 — 16H in 
the yras10 bands of medium and heavy nuclei, can be understood as a band- 
crossing between the BCS vacuum and a specific two-neutron-quasiparticle 
excited configuration that is particularly favored by the effect of rotation^ 
(see the contribution of F. Stephens and I.-Y. Lee and that of P. Ring to 
this volume) . After the understanding of this novel phenomenon, it was re- 
alized that not only the yrast band but also many excited rotational bands 
can be well described by the concept of independent quasiparticle excita- 
tions in the rotating framed This is quite nontrivial; complex rotational 
spectra at high-spin states can be nicely described by the so-called cranked 
shell modeP^ (see the contribution of S. Frauendorf to this Volume), which 
is one of the most important achievements in the studies of rapidly rotating 
nuclei. 

The Cooper pair in nucleus is composed of a pair of nucleons in the time- 
reversal conjugate orbits whose angular momenta couple to J = 0@ The 
effect of rotation, which appears as the Coriolis and centrifugal forces in 
the rotating frame, tends to align the angular momenta of nucleonic orbits 
to the rotation axis, and consequently breaks the Cooper pairs. In analogy 
to the metallic superconductors in the magnetic field, it was predicted that 
the phase transition from the supcrfluid to the normal phase is induced by 
the rapid rotation.^ However, a sharp transition as in macroscopic systems 
would not be expected in a finite system such as the nucleus. Instead, the 
finite nuclear system provides the opportunities to study a "phase transi- 
tion" in terms of the individual quantum states such as the rotational-band 
spectra with non-negligible effects of the dynamic fluctuations. In fact, 
the transition is not very simple even within the mean-field approximation: 
The effect of the band-crossings (backbendings), i.e., the successive excita- 

a The word "yrast" means dizziest, and the yrast state is the lowest energy state at a 
given angular momentum. Connecting the yrast states composes the yrast band. 
b The nucleon pairs with higher multipole, e.g. the quadrupole pair (J = 2), also play 
important roles especially in deformed nuclei. 
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tions (alignments) of quasiparticles, is more dramatic,^ and the calculated 
pairing gap reduces stepwisely along the yrast states. It is now believed 
that the unpaired phase is realized for neutrons at spins / « 20 — 30h in 
the rare-earth region, evidence for which is given by comparing the ob- 
served spectra with the rotating single-particle energies with zero pairing 
gap! 17 ! 18 ! However, it was recognized that the effects of pairing correla- 
tions remain considerably after vanishing the static (BCS) pairing gapj 19 * 20 ^ 
the "effective pairing gap" including the dynamic fluctuations beyond the 
static mean-field does not vanish and only gradually decrease across the 
phase transition.EE^ 

In the following, after briefly reviewing how to treat the nuclear ro- 
tational motion, I discuss the theoretical method to evaluate the pairing 
fluctuations within the random phase approximation (RPA) JUM] A few ex- 
amples of the calculated results, taken from our studies in Refs. Il|23rl27l 
are presented in comparison with experimental data. 

2. Description of rotational motion — Cranking model 

In order to make this article self-contained, here I recapitulate the method 
to treat the rotational motion and to analyse the rotational spectra; see 
e.g. Refs. I5I6I13I14I for detailed accounts. 

The nuclear collective motion is treated semiclassically, which is called 
the "cranking" prescription.^ Namely the Hamiltonian of the system is 
transformed into the uniformly rotating fram^l, 

H = H — LO m tJx, (1) 

where w ro t denotes the rotational frequency about the rotation axis (x-axis) , 
which is chosen to be one of the principal axes of the deformed body with 
largest moment of inertia, and is usually perpendicular to the symmetry 
axis of the quadrupole nuclear shape. Since we are mainly interested in 
the lowest energy (yrast) high-spin states, this is a natural assumption (see 
Ref. [8] for more general situations) . The energy in the rotating frame E' = 
(H') 1 which is called the "routhian", and the angular momentum along the 
rotation axis, I x — (J x ) = —dE'/doj ro t with I x = + I) « I + h, are 

evaluated as functions of the rotational frequency u; rot . 

On the other hand, the nuclear collective rotation is measured as the 
rotational spectra, E(I), which are composed of a group of states with dif- 
ferent angular momentum / changing by two units ( AI = 2), and connected 

c h = 1 unit is used for mathematical expressions. 
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by the strong electric quadrupole (E2) 7-ray emissions. In accordance with 
the simple assumption of rotational motion in Eq. ([1}, the rotational fre- 
quency is calculated by 

dE _ E(I + 1)-E(I-1) _1 

^rot(I) - ^ ~ (J+1) _ (J _ 1) - 2^< ( 2 ) 

with the 7-ray energy E y of the associated rotational transition. This 
implicitly defines the relation I x (u> IQ t), between the angular momentum 
I x and the rotational frequency w ro t, and then the experimental routhian 
E'(w to t) is obtained as 

E'(u Tot ) = E(I(u Tot )) - u mt I x (uj ro t). (3) 

In this way the theoretical routhians can be directly compared with the 
experimental routhians, although the latter are given only at the discrete 
points of the rotational frequencies. 

In the mean-field approximation, e.g., in the cranked shell model, the 
Hamiltonian H is replaced with the one-body Hamiltonian, 

H ->• h = h dcf - Z\(pt + P) - XN, (4) 

where hd e { describes the single-particle motion in the deformed average 
potential, the second term is the pair-field with P^ being the monopole 
pair creation operator, 

P^ = — c] c| (i: time reversed orbit of i), (5) 

i 

and the last term — XN ensures the correct particle number on average, 
because the number conservation is broken in the BCS treatment. By diag- 
onalizing the cranking Hamiltonian with Eq. @ the quasiparticle energies 
in the rotating frame are obtained, which can be directly compared with the 
complex rotational spectra for both even and odd nuclei}^ see the contri- 
bution of S. Frauendorf to this volume for detailed explanations. Of course, 
it can be well used with A — for the case of quenched pairing correlations, 
i.e., for the normal phase routhians. 

3. Pairing fluctuations with RPA method 

The dynamic pairing fluctuations beyond the mean-field approximation 
is induced by the two-body interaction. The simple one, the so-called 
monopole pairing force, is employed with the operator P^ defined in Eq. ([5]) ; 

H = h dcf + V, V = -~(p1p + PPi), (6) 
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with the strength G. The BCS treatment of this Hamiltonian leads to the 
one-body Hamiltonian in Eq. (j4|) with the self consistent (static) pairing gap 
A = G(P t ) m f = G(P) m f , which is nothing else but the order parameter of 
the super-to-normal phase transition. 

The fluctuations about the mean-field are calculated by diagonalizing 
the Hamiltonian © within the RPA. The induced energy gain is given by 



piRPA _ 1 
corr 2 



n a>/3 



(7) 



where u> n is the RPA eigenenergy and e a is the quasiparticle (particle or 
hole) energy in the superfluid (normal) phase. They are calculated with the 
cranking prescription as functions of o; rot to study the rapidly rotating 
nuclei. Thus, the total RPA routhian is calculated as 

E RPA = E 'ml + ^orfi E 'mf = ( h dcf ~ ^rot Jx)ml ~ G{P ] ) 2 mf . (8) 

It should be mentioned that -E^r in Eq. (0 contains the exchange energy, 
E cx = (y)mf + G(pt)iLf, w hich is found to be rather constant against 
the change of u I0 t- Note that the calculation of requires all the 

RPA eigenenergies, which amount to a few or more than ten thousands 
depending on the pairing model space. Since the convergence with respect 
to the number of solutions is slow it is important to include all of them for 
stable results, which is a numerically demanding task. A general efficient 
method to perform the calculation was developed in Ref. [T] by utilizing the 
linear response theory, and it was further improved in Ref. 1271 



3.1. Response function technique 

Generally the two-body interaction can be represented by the form of multi- 
component separable force, 

i « 

^ = -^E^A" Ql = Q P , (9) 
p=i 

with Hcrmitian one-body operators Q p and strengths x p (p = 1,2, ...,<?). 
For the monopole pairing interaction ((BJ, q — 2 and Q\ = P* + P, 1Q2 = 
pt — P, and xi = X2 = G/2. The RPA eigenvalue problem can then be 
replaced to solve the following dispersion equation, 



det 72.(a;) = 0, with 72.(o;) = [1 - P(w)x] _1 P(w), (10) 
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where the q x q matrices, TZ(u>) and R(oj), are composed of the RPA and 
the unperturbed response functions for the operators Q p , and the diagonal 
matrix % = ($ P aX P ); 



R pa {u) = 



q ;(af3)q a (a(3) q p (af3)q*(a(3) 



+ ep — u) e a + ep 



(11) 



with q p (a/3) = (a/3\Q p \0) m {. Then, by employing the adiabatic turn-on 
the interaction and the analytic property of the response function , it 
was shown that the correlation energy can be calculated by the following 
formula p3 

E%£ = —j Re[tog{det[l-iZ(iu)x]}]dw, (12) 

so that it is not necessary to explicitly solve Eq. (flQ|) . Note that the integra- 
tion is taken along the upper imaginary axis z — iuj in the complex energy 
plane, for which the integrand is a smoothly decreasing function and the 
numerical integration can be done straightforwardly. In Ref. [T] a different 
integration path is taken near the positive real axis, where the integrand is 
a oscillating function, and the numerical integration should have been done 
more carefully (see Ref. 111271 for detailed discussions). 

It is instructive to consider the following RPA pairing gapj 1 * 25 * 2 -^ 



A SFA = G /i^[(0|Pt|n)(n|P|0) + (0|P|n)<n|Pt|0}] RpA , (13) 

y n ' 

in keeping with Z\np introduced in the variation after number projection 
(NP) approach^ (gee N3.4p . However, the contribution of the zero mode, 
i.e., the symmetry recovering Nambu-Goldstone (NG) mode (the pairing 
rotatiorP), which is present in the superfluid phase, diverges because of the 
small amplitude approximation inherent in the RPA: It is natural to replace 
its (divergent) contribution to that of the mean-field; 

J[l<0|PV)| 2 + |<0|P|n)| 2 j -> (Pi) mt (P) mf = {A/Gf. (14) 

I L 1 J n=NG 

Then the difference between the squared RPA and mean- field pairing gaps, 
Z^rpa — A 2 , represents the effect of pairing vibrations, which can be calcu- 
lated by integrating the trace of the RPA response matrix Tr7£(w) without 
explicitly solving the RPA equation! 1 * 25 1 27 1 

In Fig. [II an example of the RPA correlation energy and the pairing 
gaps are shown. The mean-field pairing gap A reduces stepwisely to zero 
at the critical frequency io mt = uj c w 0.33 MeV of the super-to-normal 
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Fig. 1. The RPA and BCS mean-field (mf) pairing gaps (upper) and the RPA correla- 
tion energies (lower) for neutrons in the yrast band of 164 Er as functions of the rotational 
frequency aj ro t- The results by the variation after number projection (N P) me thod are 
included as dashed lines. Here the exchange contributions are exclude d 19 ! 25 ! both in 
^RPA an< i ^np- Taken from Ref. I27l with eliminating two irrelevant lines. 

phase transition. The first reduction at oj rot « 0.24 MeV is caused by 
the two-neutron-quasiparticle alignments (excitations) related to the back- 
bending phenomenon, where the correlation energy -E^r i s discontinuous. 
At w ro t = lo c it is continuous but its derivative, i.e., the correction to the 
alignment, 8I X = — dE^.^ /duj lot , diverges, which is a drawback of the RPA 
and one has to go beyond the RPA^ or to make smooth interpolations to 
compare with the experimental data. At these two frequencies the RPA 
gap <4rpa diverges, because one of the RPA eigenenergies goes across zero. 
It should be mentioned that E^^ is almost constant as long as the BCS 
pairing gap is sizable, while its absolute value decreases after its quench- 
ing; therefore the effect of E^.^ is important after the static pairing gap 
becomes small. In contrast to the mean-field gap, Z\rpa keeps finite values 
even at highest frequencies, reflecting that the pairing fluctuations remains 
considerably in the normal phase. These behaviors of the correlation energy 
and the pairing gaps are rather general in rapidly rotating nuclei.^ 
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3.2. Routhians and alignments in normal deformed nuclei 




h oj rot (MeV) 



Fig. 2. Calculated and experimental routhians e' (left-hand side) and alignments i 
(right-hand side) for the lowest three configurations in 168 Yb; the top panels display 
the calculation without fluctuations, middle with fluctuations, and bottom experimental 
data. Taken from Ref. [T] 

In order to discuss how the correlation energy affects rotational spec- 
tra at high-spin states, the routhians e' and the aligned angular momenta 
("alignments") i for the lowest three configurations in the nucleus 168 Yb 
are shown in comparison with experimental data in Fig. [5] Here these 
quantities are plotted relative to the so-called rigid-body reference, i.e., 



e'(Wrot) = E'(uj rot ) + - Jo W r 2 ot , i(Wrot) = 4(w ro t) - Jo Wrot, ( 15 ) 
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where J7o is the rigid-body moment of inertia. In this calculation with a 
rather simple interaction in Eq. ([6]) the experimental moment of inertia 
cannot be described correctly, and the j7o value for theoretical results is 
adjusted so as to reproduce the lower frequency part of routhians e', for 
which the correlation energy remains almost constant as it is shown in 
Fig. [TJ Apparently, the calculated routhians of the higher frequency part is 
smaller than the experimental data without the pairing fluctuations. The 
(+,0) configuration is the band with positive parity and even spins and 
corresponds to the yrast band; the kink of its routhian at w ro t ~ 0.28 
MeV corresponds to the two-neutron-quasiparticle crossing. In this band 
the mean-field pairing gap almost quenches around w rot « 0.4 MeV, and 
certainly the effect of the pairing fluctuations becomes more evident at 
larger frequencies. The negative parity excited bands with (— , 0) and (— , 1) 
are two neutron excited configurations and their pairing gaps are about 
60% of the (+, 0) band at lowest frequency. Therefore, their static pairing 
correlations are reduced more than that of the (+, 0) band, and the effects of 
pairing fluctuations are more conspicuous in the relatively lower frequency 
region. With these effects of the fluctuations, the overall agreement between 
the calculation and the experiment apparently improves. From the general 
dependence of on cj rot , the correction to the alignments i is always 

negative, which is called "dealignment" , and it amounts to 2 — 3H; again, 
this makes the agreement of alignments much better. 

3.3. Moments of inertia in superdeformed nuclei 

It was also discussed^ that the pairing fluctuations play important roles in 
the nuclei with very large deformation, which are called "superdeformation" 
and very regular rotational bands have been systematically observed; see 
Refs. 130(311 and the contribution of P.-H. Heenen to this Volume. 

For the analysis of these superdeformed bands, the two moments of 
inertia are utilized quite often; they are called the kinematic and dynamic 
inertia, and respectively, and are defined by 

m 4 1 dE' dl x d 2 E' 
J^ = — = , JW = = -— 2 . (16) 

The corrections induced by the pairing fluctuations to these inertia, r5jT < - 1 - > = 
-(l/uJiot)(dE co „/duj Iot ) and SJ^ = -d 2 E COII /d 2 u> mt , are schematically 
depicted in Fig. [3J Here the symbol uj* denotes the frequency of the inflec- 
tion point in the correlation energy and is located, in most cases, near the 
critical frequency of the pairing phase transition. 
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Fig. 3. Schematic figure depicting the (smoothed) pairing correlation energy and its 
influence on the two moments of inertia iJW and . 
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Fig. 4. Calculated and experimental moments of inertia and J^ 2 ' for the yrast 

superdeformed band in 159 Gd; thick (thin) lines denote the result with (without) pairing 
fluctuations by using the RPA method. Taken from Ref. 1241 

The superdeformed nuclei are realized by the strong shell effects based 
on the special deformations, e.g., the integer axis ratio like 2:1, and reflect 
the characteristics of the deformed closed shell. Therefore, just as in the 
case of the magic nuclei, the pairing correlations are very much reduced. Es- 
pecially, those in the mass number A « 150 region are believed to be in the 
normal phase (A — 0) already in their lowest states, and then the inflection 
point uj* in Fig. [3] is expected to be lower than the experimentally observed 
frequency region. Thus, the corrections to the inertia is negative for 
and positive for j( 2 \ In Fig. [U the calculated and experimental inertia 
are compared for the yrast superdeformed band of 149 Gd. The mean-field 
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calculation overestimates while it underestimates J^ 2 \ and a good 
agreement is obtained by including the effects of pairing fluctuations; again 
this is rather generic for superdeformed nuclei in the A « 150 regionPH 

3.4. Gauge symmetry restoration 

In the RPA the broken symmetry is signalled by the appearance of the zero- 
energy NG mode, whose contribution to the correlation energy is largest, 
ss 2 A, from Eq. ([7|). This implies the importance of restoring the gauge 
symmetry (the number conservation); the method is called the number 
projection, which explicitly projects out wave functions with good particle 
numbers from a gauge-symmetry broken wave function. In fact, it has been 
knowrJ22l that the correlations beyond the mean- field approximation can be 
taken into account by optimizing the superfluid mean-field wave function 
from which the projection is carried out; the so-called variation after pro- 
jection (VAP) approach. Therefore the variation after number projection 
(NP) is an alternative method to evaluate the pairing fluctuations in the 
rotating nuclei (see the contribution of J. L. Egido to this Volume). 

Since the expectation value of the monopole pairing operator P* van- 
ishes for the number conserving wave function, the NP pairing gap is defined 
by the following^ 

^np = G^i((Ptp) NP + (ppt) NP ). (17) 

The correlation energy and the pairing gap evaluated by the NP approach 
are also included in Fig. [T] It can be seen that both quantities behave 
quite similarly to those evaluated by the RPA method, although the NP 
correlation energy is smaller indicating that the RPA method takes more 
correlations into account. In Ref. 25 comparison of the RPA and NP meth- 
ods were performed for the routhians and alignments in rapidly rotating 
nuclei, and it was found that indeed two methods give very similar results. 
A merit of the NP method is that its result is smooth across the critical 
point of the super-to-normal phase transition in contrast to the RPA. It 
is especially useful to calculate the and J^ 1 inertias, which require 
the first and second derivatives of the correlation energy. An example is 
shown in Fig. [5j where the NP method is applied^ to the yrast superde- 
formed band in 190 Hg. The inertias of the superdeformed nuclei in the 
mass number A « 190 region systematically show increasing behaviors as 
w ro t- Because of smaller shell gaps than those in the A rs 150 regions, the 
stronger pairing correlations are expected. Namely the inflection point uj* 
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in Fig. [3] is in the higher frequency range that is experimentally observed, 
and then the increasing trends of and can be obtained by the NP 
approach without recourse to the smooth interpolation, which is necessary 
for the RPA. 




(+, 0) configuration : 7r6'W 
0.1 0.2 0.3 0.4 0.5 

ftwroi [MeV] 



Fig. 5. Calculated and experimental moments of inertia and for the yrast 

supcrdeformed band in 190 Hg by using the variation after number projection method. 
Taken from Ref. 26 but with newer experimental data. 

4. Summary 

In the this contribution I explained how the effects of the pairing fluctua- 
tions appear in rapidly rotating nuclei. By making use of the response func- 
tion technique, the correlation energy induced by the pairing fluctuations 
can be evaluated within the RPA method. The calculated RPA correlation 
energy is rather constant as long as the static pairing gap is sizable, but 
its absolute value decreases after the static gap is quenched. In this way, 
the pairing fluctuations result in dealignments of about a few units with 
respect to the mean field calculation, which makes the agreements with ex- 
perimental data much better in both normal deformed and superdeformed 
nuclei. Thus, the effects of the pairing fluctuations are important especially 
after the normal phase being realized at high-spin states. 
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